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ABSTRACT 

XWAVE  is  a computer  program  for  calculating  the  steady- 
state  pressures  in  an  infinite  fluid  exterior  to  a closed,  arbi- 
trary shaped  structural  surface  for  which  a normal  velocity 
distribution  has  been  determined.  For  this  report  the 
velocity  distribution  is  considered  to  result  from  a vibrational 
motion  of  the  surface.  The  method  of  the  program  is  based  on 
a numerical  solution  of  the  surface  Helmholtz  integral  equation 
which  relates  surface  pressure  to  velocity  and  a subsequent 
numerical  Integration  of  the  Helmholtz  exterior  field  equation 
which  relates  field  pressures  to  the  surface  velocity  and 
pressure  distribution. 

The  program  offers  a variety  of  capabilities  including 
automatic  generation  of  surface  acoustic  models  for  certain  sur- 
faces of  revolution,  automatic  generation  of  several  types  of 
velocity  boundary  condition,  an  option  for  Incorporating 
structure-fluid  Interaction  effects  through  use  of  surface 
mobility  data,  and  the  use  of  input  data  to  dynamically  allocate 
computer  core  storage  for  arrays. 

Three  calculations  Illustrate  XWAVE  data  configurations  as 
well  as  some  applications  of  the  program  to  problems  involving 
vibrating  structural  surfaces. 


ADMINISTRATIVE  INFORMATION 

This  work  was  sponsored  by  Naval  Sea  Systems  Command  (03F)  with  fund- 
ing provided  under  Subproject  SF  53  532  020,  Task  15325,  Work  Unit 
1-1808-009. 

INTRODUCTION 

This  report  documents  the  use  of  XWAVE,  a digital  computer  program 
for  computing  the  radiated  sound  pressure  from  the  surface  of  a vibrating 
structure  submerged  in  an  infinite  fluid.  The  particular  problem  con- 
sidered is  that  of  a structure  with  arbitrary  closed  surface  S whose 
points  s vibrate  with  time-harmonic  motion  as  Indicated  in  Figure  1. 

The  pressure  in  the  surrounding  fluid  then  also  varies  harmonically 
and  the  wave  equation  will  have  the  form 

V^p  + k^p  = 0 (1) 


where  p “ sound  pressure  in  the  fluid  (pressure  above  ambient) 
k “ wave  number  of  vibration;  k - ^i/c 
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Figure  1 - Vibrating  Structure  Immersed  in  an  Infinite  Fluid 

u)  = angular  frequency  of  vibration 
c “ speed  of  sound  in  the  fluid 
iw  t 

The  factor  e is  divided  out. 

The  boundary  condition  on  the  structural  surface  contacting  the 
fluid  is 

= - Iwp  V (s) 

3n  n 

where,  referring  to  Figure  1, 

s “ point  on  the  closed  surface  S of  arbitrary  shape 
n = direction  normal  to  the  structural  surface  at  s 
v^  = normal  component  of  velocity  magnitude  v at  s 
P » fluid  density 
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The  boundary  condition  at  infinity  is 

glkjzl 

p(z)  ~ ; a - » C 

where  z denotes  a point  in  the  fluid. 

XWAVE  achieves  its  solution  to  this  boundary  value  problem  by  first 
calculating  the  numerical  solution  to  a surface  Helmholtz  integral 
equation  formulation  that  relates  the  unknown  surface  pressure  dis- 
tribution to  a normal  surface  velocity  distribution  which  is  specified. 
The  field  pressures  are  solved  by  simple  quadrature  of  an  integral 
formulation  that  relates  the  field  pressures  to  the  surface  pressure 
and  normal  velocity. 

The  integral  formulations  and  associated  numerical  solution  tech- 

1* 

niques  which  the  program  uses  have  already  been  discussed.  Various 

2-4** 

practical  applications  of  the  program  have  also  been  reported. 

The  purpose  here  is  to  briefly  review  the  Integral  equations  and 
their  forms  for  numerical  solution,  to  summarize  the  types  of  data 
required  by  XWAVE,  and  to  illustrate  how  the  data  are  assembled  for 
running  applications. 
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OVERVIEW  OF  EQUATION  FORIfULATIONS  FOR  SURFACE  AND  FIELD  PRESSURE 
The  heart  of  the  XWAVE  calculation  is  the  numerical  solution  of  the 
Helmholtz  equation  for  surface  pressure: 


P(i') 

~ 


f r ® 

-^iklx'-xi' 

JJ  ^ 

4n|y'_^| 

S 


II 


9P(X) 

3n 


(4) 


where  y and  denote  surface  points. 

The  program  uses  a form  of  this  equation  in  which  both  sides  have  been 
divided  by  *^cVq  with  v^  a reference  velocity  for  the  vibration  mode  of 
the  surface  and  in  which  the  following  substitutions  have  been  made: 
Dimensionless  surface  pressure  p(^)  = p(^)/PcVq 
Dimensionless  surface  normal  velocity  v(^)  = v(y)/v^ 

Boundary  condition  from  Equation  (2) 


(5) 


The  numerical  method  of  solving  Equation  (5)  is  based  on  a subdivision 
of  the  structural  surface  S into  N interconnecting  patches  of  area  or 
stations  as  indicated  in  Figure  2.  The  kernel  functions  of  Equation  (5) 


X 


can  then  be  evaluated  for  each  station  paired,  in  turn,  with  all  stations 
including  Itself  to  yield  coefficients  of  a system  of  linear  equations 
for  the  pressure  at  each  station.  Since  details  of  the  evaluation  are 
already  available  (see  Reference  1 and  the  single  asterisk  footnote  on 
page  3),  it  is  sufficient  for  present  purposes  simply  to  recall  the  form 
of  the  equations  (see  Equation  (15)  from  CMD-2A-71): 


(lG2^,l  lA.  .1 


2 


where  J ; 

[I] 


{v.} 

J 


i.  j 

N 


[I]^  {pj  = ik  [G^j]  [Ajj]  <Vj),  1,  j = 1,  2 N (6) 

* matrices  whose  elements  are  the  evaluated  kernel 
functions  for  stations  i,j 

= identity  matrix 

= diagonal  matrix  whose  elements  are  the  station  areas 

= vector  of  nondimensional  pressures:  p.  = p,/pcv_ 

J J 0 

with  Vq  a reference  velocity  for  the  vibration  mode 

= vector  of  surface  normal  in-fluid  velocities 

nondimenslonalized;  v.  = v , /v„ 

J J 0 

refer  to  station  numbers 
= number  of  surface  stations 


The  structure-fluid  interaction  effect  can  be  Incorporated  directly 
into  Equation  (6)  through  the  modification  (see  Equation  (B5)  in 
CMD-24-71) : 


[c,jnA,,iiuj)i  i.j 


(7) 


= 1,  2,  . . . , N 


where  q^^j  is  the  in-vacuo  normal  surface  velocity  response  at  Station  i 
due  to  a unit  force  of  frequency  w acting  normal  to  the  surface  at  Station 
j,  and  is  a vector  of  the  In-vacuo  normal  surface  velocity  for  the 

vibration  mode. 
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When  the  wave  number  k is  not  at  or  near  one  of  the  set  of  Internal 

resonant  wave  numbers  of  the  cavity  enclosed  by  the  structural  surface,  a 

solution  for  surface  pressure  can  be  obtained  from  Equation  (6).  Other- 

i 3 

wise,  some  special  computing  procedure  must  be  used.  ’ A similar  con- 
dition applies  for  Equation  (7)  although  In  this  case  the  set  of  critical 
frequencies  may  differ  from  those  of  Equation  (6). 

The  surface  pressure  solution  together  with  the  normal  surface 
velocities  permit  field  pressures  to  be  calculated  from  the  integral 


formulation 


II . 


I z-\ 


5n  1 z-^ I 


For  near  or  intermediate  field  pressures,  the  matrix  equation  which  can  be 
used  to  approximate  Equation  (8)  by  using  the  subdivided  surface  S 
(Figure  2)  is  (see  Equation  (20)  of  CMD-24-71); 


(p(z^) } 


i = 1,  2, 

ik[G  ][A  ](v(j;,)}  + [G2  ] [A  ] (p(yj), 

JJ  J ■*‘J  Jj  “j  1 = 1 9 

J A » ^ » 


where  M = number  of  field  points  selected 
N = number  of  surface  stations 
= location  of  field  point 
y^  = surface  station  location 

[G^j]  = matrix  whose  elements  are  generated  by  evaluating  the  function 

y ^ ) /(4ir  I 2-^  I ) for  each  field  point  z.  paired  in  turn 
with  all  surface  stations 

[G2^j]  = matrix  whose  elements  are  generated  by  evaluating  the 
f unct  ion  ^S/(4’'|z-jf|)jfor  each  field 

point  paired,  in  turn,  with  all  surface  stations. 


A more  simplified  form  of  Equation  (8)  can  be  derived  for  the  pressure 
calculation  at  the  far-field  poi'.^s;  see  CMD-24-71.  With  reference  to  the 
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notation  In  Figure  3,  the  form  is  recalled  to  be  (see  Equation  (13)  of 
CMD-24-71) 


• SURFACE  OF  LARGE  SPHERE 
CENTERED  AT  0 


FAR-FIELD  POINT 


I 




where 


Figure  3 - Geometry  for  Far-Field  Pressure  Calculation 


(j[)  + P(^)  cos  P]e 


-ikl^lcos  a 


. 4it  I z I 

e •»’ 


The  matrix  equation  which  approximates  Equation  (10)  for  the  subdivided 
surface  S Is  (see  Equation  (23)  of  CMD-24-71): 


{F(z^)>-[Gpj,  1 [Ajj]{v(yj)}+[G2pp  HA^ ^ ) {?(^j ) } , 


j 


1.  2, 
1.  2, 


^ar-Fleld 

(11) 

N 


where 


\ar- 


Field 


% 

(G 


FF 


IJ 


[G2 


FF 


ij 


number  of  far-fleld  points  selected 
« location  of  a far-field  point 

= matrix  whose  elements  are  generated  by  evaluating  the 

function  e “ foj-  each  far-field  point 

paired  in  turn  with  each  surface  station 

= a matrix  whose  elements  are  generated  by  evaluating  the 

function  cos  6 e ® for  each  far-field  point 

paired  in  turn  with  each  surface  station  y^. 


SURFACE  GEOMETRIC  DATA 

ARBITRARY  SHAPES 

In  addition  to  the  wave  number  of  vibration  k,  various  geometric  data 
for  the  structural  surface  S,  Figure  2,  are  required  to  set  up  the  co- 
efficients of  the  matrix  equations  (Equation  (6)  or  (7))  for  surface 
pressure.  For  a particular  subdlvislonlng  of  the  surface,  one  must  specify 
the  following  data  for  each  station: 

1.  Station  number  (an  Integer  1 < !,  < N) 

2.  Surface  area 

3.  The  X-,  y-,  z-coordlnates  of  a centrally  located  point  (indicated 
by  . in  Figure  2)  within  each  station;  this  point  will  be  designated 

as  a "base"  point  for  the  station 

4.  The  X-,  y-,  z-coordlnates  of  unit  normal  vector  to  the  surface  at 
the  base  point 

5.  Local  curvature  of  the  surface  with  reference  to  the  base  point 
If  the  fluid-surface  Interaction  effect  is  to  be  calculated  by 

using  Equation  (7),  the  above  data  are  extended  to  Include  the  sound 
speed  in  the  fluid  c,  the  fluid  mass  density  p,  and  the  elements  q^j  of 
the  surface  dynamic  mobility  matrix  (see  page  5).  The  q^j  are  assumed 
here  to  be  available  from  a source  external  to  XWAVE  and  may  derive  from 
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I 


ft 


an  analytic,  numerical,  or  experimental  analysis  of  surface  mobility.  The  I 

program  requires  only  that  these  data  be  compatible  with  the  particular  J 

surface  modeling  being  used  and  that  they  are  in  the  proper  format  as  ? 

/ 

described  In  the  section  entitled  DATA  FORMATS.  i 

XWAVE  incorporates  several  options  for  supplying  the  data  indicated  I 

in  Items  2-5  to  the  subroutines  which  calculate  matrix  coefficients. 

One  option,  the  default  option,  is  used  with  those  surfaces  for  which 
automatic  generation  of  the  above  data  has  not  yet  been  included  in  the 
program.  For  such  cases  the  entire  geometric  data  for  the  surface  model 
are  specified  on  punched  cards  to  be  read  into  storage  by  the  program. 

The  remaining  options  for  surface  geometric  data  pertain  to  certain 
surfaces  of  revolution  which  have  been  used  in  testing  and  applying  the  ' 

program. 

SURFACES  OF  REVOLUTION 

XWAVE  presently  has  the  capability  for  automatically  generating  the 
surface  geometric  data  for  three  types  of  revolved  surface.  We  shall 
consider  in  detail  the  data  specified  to  the  program  for  each  surface.  ? 

Piecewise  Conical  Shell 

These  are  defined  as  surfaces  formed  by  revolving  curves  con-  ! 

slsting  of  straight  line  segments  joined  end  to  end.  Examples  of  this  ; 

kind  of  surface  are  shown  in  Figure  4.  j 

The  primary  portion  of  surface  utilized  in  generating  the  geometric 
data  is  defined  here  to  be  a "region."  A region  is  formed  by  revolving  J 

either  the  full  or  partial  length  of  a straight  line  from  the  generating 

curve  through  a specified  arc.  The  data  required  to  specify  the  extent 
of  a region  to  the  program  are  indicated  in  Figure  5.  In  the  figure, 

Zj,  z^  indicate  the  extent  of  the  region  along  the  z-axis,  r^  deter- 
mine the  line  segment  slope,  and  6^,  6^  indicate  the  extent  of  the 

region  in  6,  6^  and  62  being  measured  from  the  positive  y-axls  j 

counterclockwise.  ; 

J 
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Figure  A - Surfaces  Generated  by  Revolving  Curves  Formed  f 

Straight  Line  Segments 


Figure  5 - Geometry  of  Surface  Region 


formed.  The  disk  represents  the  case 

Zj  . Z2 

. tj  - 0,  r^  > 0 

= 0 degrees,  9^  = 360  degrees 


O G (p 

DISK  PARTIAL  DISK  CYLINDER  CYLINDRICAL 

SECTION 


ANNULAR  RING  PARTIAL  ANNULAR  CONE  PARTIAL  CONE 

RING 


RIGHT  CONICAL  PARTIAL  SECTION  OF 

SECTION  RIGHT  CONE 

Figure  6 - Typical  Surface  Regions  from  Revolving  Straight  Line  Segments 
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and  the  annular  ring  represents  the  case 


> 0,  r2  > r^ 


* 0 degrees,  6^  • 360  degrees 


The  subdivlsionlng  of  the  surface  within  a region  into  stations  Is 
controlled  by  two  Integers  n and  m supplied  as  data.  The  first  integer 
designates  the  number  of  equal  increments  (Az  - {z^-z^)(n)  into  which 
the  region  interval  in  z is  to  be  divided.  When  projected  to  the  surface, 
these  intervals  define  n bands  of  area  as  indicated  in  Figure  7. 


Figure  7 - Surface  Region  Subdivided  into  Bands  of  Area 


i 


The  second  Integer  m designates  the  number  of  equal  Increments  into 
which  the  6 interval  of  the  region  is  to  be  divided;  A6  » 

The  pair  n,m  then  controls  the  station  modeling  of  the  region  as  shown 
in  Figure  8. 


X 
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Figure  8 - Surface  Region  Subdivided  into  Stations 

The  XWAVE  data  generator  automatically  locates  a "base"  point  within 
each  station  area  at  a position  1/2  AZ,  1/2  A6.  Coordinates  of  a unit 
normal  and  the  mean  surface  curvature  are  calculated  for  the  "base" 
point  by  using  formulas  previously  discussed  (see  single  asterisk  footnote 
on  page  3).  If  the  generating  line  for  the  region  is  perpendicular  to  z, 
the  coordinates  of  the  unit  normals  to  base  points  within  are  simply  x=0, 
y=0,  z=  +1.  The  sign  of  z is  resolved  by  specifying  it  as  an  item 
of  data  accompanying  the  limits  on  extent  of  the  region  (page  9)  and  m 
and  n.  The  sign  is  easily  resolved  for  all  other  regions  from  the  slope 
of  the  generating  line. 

The  surface  curvature  used  at  each  station  is  taken  to  be  the  mean 
curvature  at  the  base  point.  The  mean  curvature  is  defined  as  one-half 
the  sum  of  the  maximum  and  minimum  principal  normal  curvatures  at  the 
base  point.  For  surfaces  of  revolution,  these  principal  lines  of 
curvature  are  the  meridian  and  parallel  passing  through  the  base  point. 
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Note  from  Figure  6 that  the  surface  regions  under  consideration  are 
either  flat  or  have  curvature  only  In  the  direction  of  the  parallels.  In 
the  first  Instance,  the  curvature  Is  obviously  zero  for  all  stations 
within  the  region.  In  the  second  Instance,  curvature  at  the  base  point 
Is  the  curvature  of  the  local  parallel  which  Is  given  by  the  distance 
from  the  base  point  to  the  axis  of  revolution  along  a line  perpendicular 
tc  the  surface. 

The  area  of  each  station  Is  calculated  as  I/mth  the  area  of  the  band 
In  which  It  Is  located.  The  numbering  of  stations  Is  according  to  the 
order  In  which  they  are  generated  by  the  program;  this  is,  counterclockwise 
about  each  band  and  running  successively  through  all  bands  from  the  left 
to  right  end  of  the  structural  surface  as  oriented  In  Figure  A.  Figure  9 
Illustrates  the  ordering  of  stations  for  three  adjacent  bands  of  area  In 
two  possible  modeling  situations.  The  second  of  these  (b)  represents 
the  case  where  a surface  region  360  degrees  in  extent  (for  example  a 


X X 


X 


Figure  9 - Ordering  of  Surface  Station  Numbers  for  Bodies  of  Revolution 
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disk  or  cone  as  in  Figure  6)  and  containing  a single  station  is  used  in  i 

closing  off  the  end  of  a structural  surface.  The  dots  in  Figure  9 1 

designate  station  base  points.  As  indicated,  the  reference  station  for  j 

t 

beginning  the  numbering  in  each  band  (except  end  bands  with  one  station) 
is  the  first  station  generated  above  the  positive  y-axis. 

Prolate  Spheroid  j 

Here  surfaces  are  formed  by  revolving  ellipses  with  major  axis  along 
the  axis  of  revolution.  The  surface  profile  is  specified  to  the  data  ■ 

I 

generator  through  two  input  parameters  AA  and  BB;  these  are  respectively  i 

the  major  and  minor  axes  of  the  generating  ellipse  which  is  assumed  to  be  j 

centered  about  the  origin  of  the  coordinates.  i 

Figure  10  shows  a typical  region  on  this  surface.  i 


i 


POLE 


Figure  10  - Typical  Region  of  a Prolate  Spheroidal  Surface 


To  specify  a prolate  spheroidal  region,  the  following  data  are 
supplied  to  the  data  generator: 

1.  6p  the  angular  extent  of  the  region  with  respect  to  the 

generating  axis  z;  6 Is  measured  counterclockwise 
from  the  positive  y-axls;  0 degrees  < 6 < 

360  degrees 

2.  angular  extent  of  the  region  In  colatltute  measured 
from  the  z-negatlve  pole  of  the  body;  0 degrees 

< ♦ < 180  degrees 

Modeling  of  the  surface  region  Is  determined  by  Integer  data  n and  m 
for  the  region.  The  first  Integer  specifies  the  subdlvlslonlng  of  the 
region  Into  n equal  bands  of  latitude  (A^  « (4'2-0j)/n  degrees);  the  second 
Integer  specifies  m equal  strips  of  longitude  (A6  = (02-9^)/m  degrees), 
as  Indicated  In  Figure  11. 

The  area  of  each  station  within  a band  of  latitude  A0  is  1/mth  total 
area  of  the  band.  The  base  point  within  each  band  is  located  at  1/2  A6, 
1/2  Alt. 


Figure  11  - Surface  Modeling  of  Prolate  Spheroidal  Region 


A mean  curvature  for  each  station  is  found,  as  in  the  case  of  coni- 
cal shell  regions,  from  the  curvatures  of  the  meridian  and  parallel  through 
the  base  point;  the  former  curvature  is  readily  determined  from  the 
equation  for  the  particular  ellipse  in  the  xz-plane. 

The  implied  numbering  of  stations  generated  follows  the  pattern 
illustrated  in  Figure  9;  single  end  stations  as  indicated  in  part  (b)  of 
the  sketch  now  become  spheroidal  caps. 

Sphere 

The  particular  spherical  profile  is  specified  to  the  data  generator 
through  a single  input  parameter  RADIUS.  The  data  to  specify  surface 
regions  and  their  subdivision  into  stations  are  the  same  as  just  described 
for  the  prolate  spheroid. 


SYMMETRY 


It  can  easily  be  shown  that  when  symmetry  planes  occur  simultaneously 
for  both  the  surface  geometry  and  velocity  boundary  conditions,  they  are 
also  symmetry  planes  for  the  surface  pressure  solution.  In  other  words, 
if  the  surface  has  a symmetry  plane  which  is  also  a symmetry  or  anti- 
symmetry plane  for  the  boundary  condition,  then  the  plane  is  respectively 
a symmetry  or  antisymmetry  plane  for  the  pressure  solution.  This  result 
enables  one  to  specify  some  basic  portion  of  the  total  surface  for 
modeling  and  to  obtain  the  modeling  of  the  remainder  of  the  surface 
through  reflections  of  the  basic  model  through  symmetry  planes. 

XWAVE  incorporates  options  for  using  the  principal  coordinate  planes 
as  symmetry  planes  in  specifying  the  structural  surface  model.  The 
elements  or  stations  which  model  the  nonredundant  portion  of  the  sur- 
face with  respect  to  these  symmetry  planes  are  referred  to  as  "basic" 
elements.  The  elements  which  model  the  remainder  of  the  surface  are  ob- 
tained by  reflecting  the  basic  elements  in  the  symmetry  planes.  Figure  12 
Illustrates  a basic  surface  element  and  the  images  of  successive  reflections 
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The  total  of  basic  plus  reflected  elements  gives  the  "effective"  | 

number  of  elements  for  the  surface  model.  The  total  computer  core  S 

storage  required  for  the  surface  geometric  data  is,  however,  only  that  i 

required  for  the  basic  elements  since  the  data  for  reflected  elements  | 

can  be  obtained  from  the  basic  element  data  by  appropriate  sign  changes. 

With  regard  to  generation  of  the  coefficient  matrix  (for  example. 

Equation  (6))  for  surface  pressure,  the  program  recognizes  principal 
coordinate  plane  and  axial  symmetry.  In  the  case  of  principal  coordinate 


plane  symmetry  alone,  the  order  of  the  matrix  equation  is  equal  to  the 
number  of  basic  elements  or  stations  used  to  model  the  nonredundant 
portion  of  the  structural  surface.  The  ijth  coefficient  of  the  matrix 
is  composed  of  a sum  of  terms,  the  first  of  which  represents  an  eval- 
uation of  the  kernel  function  for  the  pairing  of  basic  element  i and 
basic  element  j.  The  remaining  terms  of  the  sum,  respectively,  represent 
an  evaluation  of  the  kernel  for  pairings  of  the  same  basic  element  1 
with  all  reflected  images  of  the  same  basic  element  j in  the  symmetry 
planes.  The  sign  of  a term  is  plus  if  the  image  element  is  from 
reflection  across  a symmetry  plane  of  the  boundary  condition;  the  sign 
is  minus  if  across  a plane  of  antisymmetry. 


? 


When  the  surface  geometry  and  boundary  conditions  are  axially 
symmetric,  the  pressure  solution  need  be  obtained  for  only  one  station 
in  each  latitude  band.  We  arbitrarily  select  these  stations  on  the  x- 
posltlve  side  of  the  yz-plane  and  bordering  e * 0 degrees  (see  Figure  8). 

The  matrix  coefficients  are  generated  by  a modification  of  the  pro- 
cedure described  above  for  principal  coordinate  plane  symmetry.  In  the 
first  part  of  the  calculation,  the  index  1 is  constrained  to  reference 
just  the  single  selected  station  in  each  band  of  latitude  and  the  kernel 
function  is  evaluated  for  all  station  pairs  i,j;  reflected  images  of  J in 
the  coordinate  planes  are  accounted  for  as  before.  In  the  second  and  last 
stage,  those  terms  resulting  from  the  pairing  of  station  i with  circum- 
ferential stations  j within  individual  bands  are  added  together.  This 
addition  accounts  for  the  fact  that  the  circumferential  stations  j (where 
j 4 i)  are  actually  images  of  the  base  station  j = 1,  when  it  is  success- 
ively reflected  in  symmetry  planes  inclined  at  appropriate  angles  to  the 
coordinate  planes.  The  matrix  of  coefficients  generated  in  this  way  is 
seen  to  be  of  order  equal  to  the  total  number  of  latitude  bands  prescribed 
for  the  surface  model. 

The  procedures  described  for  generating  the  coefficient  matrix  for 
cases  involving  symmetry  apply  equally,  of  course,  to  the  generation 
of  elements  of  the  matrix  product  [G^^](Aj^]  in  the  surface  pressure 


equat ion . 


POUNDAPY  CONDITION  DATA 

The  boundary  condition  data  which  must  be  supplied  to  the  XTIAVE  pro- 


gram consist  of  the  elements  of  the  vector  (Vj } (see  Equation  (6))  or 
the  vector  {u.}  (see  Equation  (7)).  For  the  latter  case  in  which 
structure-fluid  interaction  effects  are  being  incorporated  into  the 


calculation,  the  vector  is  assumed  to  be  generated  external  to 


XWAVE,  most  likely  by  a program  for  dynamic  structural  analysis.  The 


elements  of  1^^}  are  assumed  to  be  stored,  along  with  the  mobility 


matrix  elements  on  an  auxiliary  storage  device  (disk)  or  in  card 


form  for  reading  by  XWAVE. 
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XWAVE  Incorporates  several  options  for  specifying  the  in-fluid  velocity 
vector  fVj}.  One  option  allows  the  user  to  specify  a single  value  of 

velocity  V to  be  applied  at  each  station  within  a region  on  a surface 
of  revolution  (see  Figure  8).  Assignment  of  v to  each  station  Is 
handled  automatically  by  the  program.  In  the  case  of  axial  symmetry,  v 
is  assigned  to  ot\ly  one  station  In  each  band  of  latitude  within  the  region 
(refer  to  page  19,  top). 

A second  option,  also  designed  for  use  with  surfaces  of  revolution, 
automatically  generates  a normal  velocity  distribution  from  the  expression^ 


V = v„ii)(z)  cos  Me  (12) 

n u 

where  v^  - constant  velocity  amplitude 

(j/(z)  = arbitrary  dimensionless  function  of  z 
M = integer  constant 

A variety  of  assumed  velocity  distributions  corresponding  to  vi- 
bration modes  of  surfaces  of  revolution  can  be  obtained  from  the  above  ex- 
pression. For  example,  longitudinal  modes  and  circumferentially  symmetric 
modes  with  M = 0,  flexural  modes  with  M = 1,  and  loba’-  modes  with  M > 1.^ 
Additional  distributions  can  be  obtained  from  lint-'ar  combinations  of 
these  modes. 

The  data  supplied  to  XWAVE  with  this  option  are  the  constants  v^  and 
M and  for  each  surface  region  the  values  i(z^),  i = 1,  2,  3,  ...,  LB; 

LB  denotes  the  number  of  latitude  bands  Into  which  the  region  is  subdivided. 
The  points  z^  coincide  with  the  z-coordlnates  of  the  base  point  locations 
(i.e.,  1/2  Az;  see  Figure  8)  within  each  band.  Similarly,  the  particular 
values  of  e,  9^  used  in  computing  v^  correspond  to  the  angular  locations 
(1/2  A6)  of  base  points  within  each  band. 

A third  option  allows  the  user  to  explicitly  specify  the  normal 
velocity  to  be  assigned  each  station  in  the  surface  model.  This  option 
provides  for  cases  where  the  facilities  of  the  first  and  second  option  do 
not  apply.  The  velocity  data  with  this  option  are  read  into  the  computer 
from  punched  cards. 

^Chertock,  C. , "Sound  Radiation  from  Vibrating  Surfaces,"  Journal  of 
the  Acoustical  Society  of  America,  Vol.  36,  No.  7,  pp.  1305-1313  (Jiil  196A). 
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DATA  FOR  ITERATIVE  SOLUTION  OF  THE  SURFACE  PRESSURE  EQUATION 
If  In  Equation  (6)  ve  denote 


(lC2,jll»jjl  -i  m)  i [M|,l 


and 


the  equation  can  be  written 


IV 


'"ij>  * I'lj’  '^i'  ■ 'L> 


(I  3) 


A similar  expression  can  be  written  for  Equation  (7).  XWAVE  uses  a 
Cuass-Seidel  iteration  scheme  for  solving  Equation  (13).  The  form  of  the 
iteration  is 


(n+  1 ) 


( n) 


Pi 


= P.- 


+ H . 


1 

M.  . 

1 1 


F,  - 
1 


i-1 
.1  = 1 


(n+1) 


M.  . p. 
J-J  J 


N 

y V 

j=i 


(n) 


1.1 


with  H the  relaxation  factor. 


The  initial  value  of  p^  is  given  by 

(0) 

Pi  = 


Ni 


(U) 


(15) 


The  iteration  is  terminated  either  when 


n+ 1 

iPi 


(16) 


or  when  a prespecified  bound  on  the  number  of  iterations  is  reached. 

If  the  upper  bound  on  iterations  is  met  first,  the  last  vector  computed  is 

n+ 1 n 

printed  out.  In  this  case  the  magnitude  (p 


p, I which  is  printed 
' i ' max 


after  each  iteration  can  be  used  to  estimate  the  degree  to  which  the 
the  pressure  vector  has  converged. 
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Data  Input  for  the  iterative  solution  will  Include  the  relaxation 
factor  H,  the  convergence  criterion  e,  and  an  upper  bound  on  the  number 
of  iterations. 

FIELD  POINT  DATA 

XWAVE  Incorporates  an  option  for  either  direct  specification  or 
automatic  generation  of  field  point  locations  (see  page  6)  and  ^ 

(see  page  8).  As  noted  earlier,  far-fleld  points  ^ are  assumed  to  lie 
on  a sphere  of  large  radius  centered  on  the  submerged  structure.  Near 

and  Intermediate  points  z,  lie  between  the  structural  surface  and  those 

■*►1 

distances  which,  in  an  acoustic  sense,  can  be  considered  to  be  far  field. 

To  directly  fpecify  far-field  points,  one  supplies  on  punched  cards 
the  longitude  and  latitude  of  each  point  desired  according  to  the  format 
given  in  the  section  entitled  DATA  FORMATS.  The  coordinate  system  used 
in  referencing  longitude  ♦ and  colatitude  is  shown  in  Figure  13. 


Figure  13  - Spherical  Surface  Coordinate  System  Centered  in  a Structure 

Alternatively,  one  can  select  an  option  to  automatically  generate 
points  (<!',&)  on  the  spherical  surface  at  infinity.  When  this  option  is 
used,  the  data  specified  Include: 
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1.  The  Uniting  nunbor  of  degrees  of  colatitude  extending  from 
Ip  = 0 degrees  (0  degrees  < i < 180  degrees) 

2.  The  limiting  number  of  degrees  of  longitude  extending  from 
<1  = 0 degrees  (0  degrees  < < < 360  degrees) 

3.  The  number  of  colatitudes  desired 

4.  The  number  of  longitudes  desired  vrlth  each  colatitude  other 
than  0 or  180  degrees 

An  arbitrary  distribution  of  near  or  internediate  field  points  is 
obtained  by  specifying  the  x-,  y-,  z-coord inates  of  each  point.  For 
automatic  generation  of  field  points  on  a spherical  surface,  one  specifies 
in  place  of  Items  1 and  2 above,  the  colatitude  limits  the 

longitude  limits  4>2»  data  of  Items  3 and  4,  and  radius  R 

of  the  sphere. 

DATA  FOR  OYNAMIC  CORE  ALLOCATION 

The  remaining  data  required  for  XWAVE  specify  the  sizes  of  certain 
arrays  used  by  the  program.  These  data  have  two  purposes:  (1)  to  enable 

the  size  of  arrays  which  primarily  store  input  data  (such  as  the  surface 
geometric  data,  field  points,  boundary  conditions,  etc.)  to  be  tailored 
to  meet  Just  the  requirements  of  the  particular  application  and  (2)  to 
allow  the  user  to  decide  how  much  core  is  to  be  committed  for  matrix 
manipulations,  or  in  other  words,  to  what  degree  the  matrix  manipulations 
will  be  handled  as  in-core  or  out-of-core  operations.  The  flexibility 
offered  in  transforming  from  essentially  in-core  to  out-of-core  con- 
siderably increases  the  size  of  a Job  (l.e.,  with  respect  to  number  of 
surface  elements)  that  can  be  treated.  The  ability  to  reduce  Internal 
core  requirements  also  offers  the  user  an  advantage  for  applications 
Involving  relatively  small  running  times;  by  virtue  of  the  reduce'  core, 
he  can  select  a higher  running  priority  and  better  throughput  time. 

The  economics  of  running  the  program  are  also  enhanced  by  the  core  ad- 
justment feature. 
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XWAVE  arrays,  that  la,  those  of  significant  length,  are  controlled 
through  the  use  of  14  dimensions.  The  rules  for  assigning  values  to  these 
dimensions  are  tabulaten  oclow. 


Dimension 

D1 


D2 


Value 

Total  number  of  surface 
elements  in  the  surface 
model . 


Total  number  of  surface 
elements  used  to  model 
the  portion  of  structural 
surface  which  is  nonre- 
dundant  with  respect  to 
principal  coordinate 
planes  of  symmetry. 

Note:  For  surfaces  of 

revolution  modeled  by 
using  the  automatic  in- 
put facility,  this  di- 
mension is  equal  to  the 
total  sum  of  n X m (see 
pages  12,13)  over  all 
surface  regions. 

Total  number  of  surface 
elements  in  that  portion 
of  the  surface  which 
is  nonredundant  with  re- 
spect to  axial  symmetry. 
This  dimension  is  equal 
to  the  sum  of  n (see  pages 
12,13)  over  all  sur- 
face regions. 

Maximum  number  of  rows  of 
each  matrix  ] and 

(see  page  21)  that  will  be 
stored  in  core  at  any  one 
time . 

Note:  D2  = the  order  of 

the  matrix;  corresponds  to 
totally  in-core  matrix 
operations . 


Application  Condition 

Velocity  boundary  condition 
has  no  principal  coordinate 
plane  symmetries. 

Structural  surface  geometry 
has  no  principal  coordinate 
plane  symmetries. 

Structural  surface  geometry 
and  velocity  boundary  con- 
dition have  the  same 
principal  coordinate  plane 
symmetries  but  are  not 
axially  symmetric. 


Structural  surface  geometry 
and  velocity  boundary  con- 
dition are  axially  sy- 
mmetric . 


All 
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Value 


i 


D Imenslon 
D3 


D4 


D5 


f)6 


n? 


Total  number  of 
elements  in  the 
model . 


surface 

surface 


Total  number  of  surface 
elements  used  to  model  the 
portion  of  structural 
surface  which  is  nonrc- 
dundant  with  respect  to 
principal  coordinate 
planes  of  symmetry. 

Total  number  of  surface 
regions  (see  page  9) 
used  to  model  that 
portion  of  structural 
surface  which  is  nonre- 
dundant  with  respect  to 
principal  coordinate 
planes  of  symmetry. 

Number  of  radial  points 
of  longitude  (see  page 
23,  Item  4)  defining 
selected  field  points. 

Maximum  number  m (see 
page  13)  specified  for 
any  region  of  the  sur- 
face model . 


Maximum  of  the  two 
values  defined  above. 


Same  as  D5 


Maximum  value 
(see  page  12) 
specified  for 
region  in  the 
model . 


of  n 

which  is 
any 

surface 
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Application  Condition 

Neither  the  velocity  bound- 
ary condition  nor  the 
structural  surface  geometry 
have  principal  coordinate 
plane  symmetry. 

Structural  surface  geometry 
and  velocity  boundary  con- 
dition have  the  same 
principal  coordinate  plane 
symmetries  and  may  or  may 
not  be  axially  symimetric. 

Structural  surface  geometry 
is  axially  symmetric. 


Structural  surface  is  not 
axially  symmetric  and  the 
near-field  pressure  calcu- 
lation has  been  selected. 

Structural  surface  is 
axially  symmetric  and  only 
surface  pressures  or  sur- 
face pressures  and  far- 
field  pressures  are  being 
calculated. 

The  structural  surface  _1^ 
axially  symmetric  and  sur- 
face and  near-field  pres- 
sures are  being  calculated. 

Structural  surface  geometry 
is  axially  symmetric  and  the 
velocity  boundary  condition 
given  by  Equation  (12) 

(page  20)  has  been  selected. 

Same  as  the  above. 


f 


Dimension 


Value 


)lication  Condition 


Same  as  D1 


Surface  dynamic  mobility- 
data  (see  page  5) 

are  being  incorporated  in 
the  surface  pressure  calcu- 
lation. 

Same  as  the  above. 


Maximum  number  of 
columns  of  the  matrix 
product 

Equation  (7))  that  will 
be  stored  in  core  at  any 
one  time. 

Note:  D9  = the  order  of 

the  matrix;  corresponds 
to  totally  in-core  matrix 
operations . 

Same  as  D2 

Same  as  D3 

Total  number  of  near- 
field points  selected. 

Number  of  colatitudes 
(see  page  23,  Item  3) 
selected  in  defining 
far-field  points. 

Number  of  longitudes 
(see  page  23,  Item  4) 
selected  in  defining 
far-field  points. 


The  user  specifies  only  those  dimension  parameters  appropriate  to  a 
particular  application.  In  other  words,  when  no  one  of  the  criteria  for 
a dimension  parameter  is  met,  it  is  considered  blank  or  zero  in  the  data 
input.  The  parameters  will  always  include  as  a minimum  the  set  Dl,  D2, 
and  D3.  Additional  parameters  will  depend  on  the  type  of  surface  geometry 
involved,  velocity  boundary  data,  and  whether  field  pressures  are  calcu- 
lated. For  example,  if  surface,  near-field,  and  far-field  pressures 
are  to  be  computed  using  an  assumed,  In-fluid  velocity  boundary  con- 
dition with  a nonaxlally  symmetric  surface,  the  parameters  would  be  Dl, 

D2,  D3,  D5,  D12,  D13,  and  D14. 


Same  as  the  above. 

Same  as  the  above. 

Near-field  pressure  calcu- 
lation has  been  selected. 

Far-field  pressure  calcu- 
lation has  been  selected. 


Same  as  the  above. 


Data  formats  for  specifying  dimension  parameters  are  given  In  the 
section  entitled  DATA  FOEIMATS. 

OPTION  FOR  INDUCED  MASS  CALCULATION 
Once  the  surface  pressure  distribution  corresponding  to  a 

particular  velocity  boundary  condition  {v^}  has  been  determined, 
various  quantities  of  interest  can  be  calculated  apart  from  field 
pressures.  One  of  these  quantities  is  the  induced  mass  for  a prescribed 
vibrational  mode  of  a body. 

The  body  is  considered  to  be  surrounded  by  a fluid  which  tends  to 
infinity,  there  to  be  limited  by  a closed,  rigid  spherical  surface.  The 
basis  for  computing  induced  mass  is  an  expression^  for  total  energy 
increment  W in  the  fluid  as  the  vibration  velocity  goes  from  zero  to  its 
maximum  value 


II 

s 

where  S = closed  surface  of  the  body 
V = normal  surface  velocity 
t(p)  = imaginary  part  of  surface  pressure 
u)  = angular  frequency  of  vibration 


(17) 


If  the  fluid  is  assumed  incompressible  (c  " and  k -►  0),  the  total 
energy  is  entirely  kinetic 


W = T 


t,  (p  ) V dS 


(18) 


and  equals  the  maximum  kinetic  energy  of  the  fluid. 

Substituting  nondimenslonal  forms  for  pressure  and  velocity, 
V (page  4),  into  Equation  (18),  one  obtains 


T 


2(ij 


rr  - - 1 2 

r , rr  - - 1 

JJ  Mp)v  dS  - 2 Vq 

r* 

p/k  JJ  l(p)v  dS 

^ C 

p and 


(19) 


The  bracketed  quantity  is  Interpreted  as  the  Induced  mass  of  the  vibration 
mode. 
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When  the  option  for  Induced  mass  is  invoked,  XWAVE  provides  a 
numerical  approximation  to  the  bracketed  quantity  in  Equation  (19): 


M , A 

V 


Induced  mass 


/ ^ 

(p/k  i(Pi>  . . AS^ 

\ i-1 


= number  of  basic  surface  elements  with  respect  to  principal 
coordinate  symmetry  or  antisymmetry  planes 

= fluid  density 

* integer  multiplier  dependent  on  the  number  of  principal 
coordinate  symmetry  planes 

“ area  of  the  ith  surface  element 


A detailed  test  of  the  Induced  mass  capability  of  XWAVE  for  the  case 

A 

of  rigid  body  vibration  of  cylinders  has  been  documented.  A typical 
example  from  that  application  is  presented  in  the  section  entitled 
EXAMPLE  CALCULATIONS. 


OPTION  FOP  SELF-RADIATION  IMPEDANCE  CALCULATION 
A second  quantity  of  Interest  that  can  be  calculated  by  using  the 
surface  pressure  distribution  is  the  self-radiation  impedance  of  a vi- 
brating surface.  The  vibrating  surface  may  be  the  entire  closed  sur- 
face of  a body  or  a portion  of  it,  the  remainder  being  considered  rigid. 

The  self-radiation  Impedance,  which  will  be  denoted  by  is 

the  acoustic  radiation  force  per  unit  velocity  acting  on  the  portion  of 
surface  that  is  vibrating 


F /v 

ir  1 


* 

See  the  double  asterisk  footnote  on  page  3. 
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where 


dSj  (22) 

and 

“ vibrating  surface 

Vj  = surface  normal  velocity  of  the  vibrating  surface 
p,  = complex-valued  pressure  on 

The  normalized  form  of  Equation  (21)  used  for  the  calculation  is 


Z 


11 


= F / dcA  V 

ir  1 1 


(23) 


where  is  the  area  of  the  vibrating  surface  Substituting  the 

nondlmenslonal  pressure  (p^  = p^/pcv^)  into  the  expression  for 
then  yields 


Z 


11 


dSj/Aj 


(24) 


When  the  option  for  radiation  impedance  is  designated,  XWAVE  provides 
the  following  numerical  evaluation  of  Equation  (24): 


/ ^ 
' i=l 


(25) 


where  N is  the  number  of  basic  surface  elements  (with  respect  to  principal 
coordinate  symmetry  or  antisymmetry  planes)  representing  the  vibrating  sur- 
face Sj  and  ASj  is  the  area  of  the  ith  surface  element  of  Sj. 

This  option  of  XWAVE  has  been  applied  to  the  calculation  of  the  self- 
radiation  impedance  of  a ring-shaped  piston  vibrating  in  a spheroidal 
baffle.  An  example  from  this  series  of  calculations  is  given  in  the 
section  entitled  EXAMPLE  CALCULATIONS. 
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DATA  FORMATS 

This  section  summarizes  all  XWAVE  data  relating  to  calculations 
Involving  vibrating  structural  surfaces  and  Indicates  how  the  data  are 
prepared  for  reading  by  the  program.  Data  Items  which  have  been  dis- 
cussed In  detail  In  prior  sections  of  this  report  are  referenced  to  the 
appropriate  pages. 

Data  card  type  designations  Indicated  by  " " refer  to  the  data 

input  forms  (see  the  appendix)  which  have  been  prepared  for  use  with 
XWAVE. 


• "PROBLEM  TITLE"  CARD 

Columns  Contents 

1-12  Blank 

13-72  General  title  associated  with  the  run 

• "CASE  TITLE"  CARD 

Columns  Contents 

1-12  Blank 

13-72  Subtitle  associated  with  the  run 


M 


DIMENSIONS" 
1-  4 
5-  8 
9-12 
13-16 
17-20 
21-24 
25-28 
29-32 
33-36 
37-40 
41-44 
45-48 
49-52 
53-56 


CARD  (pages 
D1 
D2 
D3 
D4 
D5 
D6 
D7 
D8 
D9 
DIO 
Dll 
D12 
D13 
D14 


24-26) 


30 


3 

) 


"MISCELLANEOUS  DATA"  CARD 


I 


Columns 

Contents 

Description 

• 

1-  8 

k 

Wave  number  of  the  sound  pressure 
wave,  •=  u)/c 

9-16 

Blank 

17-24 

c 

Speed  of  sound  in  the  fluid 

25-32 

P 

Fluid  density 

33-40 

R (H) 

Real  part  of  relaxation  factor  H 

41-48 

«-  (H) 

Imaginary  part  of  H 

i 

49-5  2 

iteration  limit 

Upper  bound  on  the  number  of 
iterations  to  be  carried  out 
in  the  solution  for  surface 
pressures  (page  21) 

• 

53-60 

convergence  criterion 
= e 

A number  through  which  the  user 
may  select  the  precision  of 

the  surface  pressure  solution, 
assuming  of  course  that  a con- 
vergent solution  exists  for 
the  particular  application 
(page  21) 


• "PROGRAM  OPTIONS"  CARD 


Columns 
1-  A 


Contents 


0000 


0001 


5-  8 0000 

0001 


Description 

OPT  1 : Selector  for  structure- 

fluid  interaction  effect 

Effect  is  not  to  be  included  in 
the  calculation  and  the  normal 
surface  in-fluid  velocity  is  to 
be  specified 

Effect  is  to  be  included  and 
the  surface  in-vacuo  mobility 
coefficients  and  in-vacuo  sur- 
face normal  velocity  for  the 
vibration  mode  are  input 
(pages  5,8,9,19) 

0PT2:  Structural  surface  is 

not  a surface  of  revolution 

Structural  surface  is  a surface 
of  revolution 


31 


1 


9-12 


0000 

0001 

0002 

0003 

13-16 


0000 

0001 

0002 

0003 

0004 

17-20 


OPT3:  Selector  for  surfaces  of 

revolution 

Structural  surface  is: 

Piecewise  conical  shell  (pages 
9.10) 

Cylinder 

Sphere 

Ellipse 

0PT4:  Selector  for  surface 

velocity  distribution 
Velocity  distribution  is: 

Arbitrary  and  will  be  read  in 
from  cards  for  each  station 

Uniform  over  entire  surface 

Of  the  form  Vq'I'(z)  cos  M0 

(page  20) 

Generated  for  rigid  body  motion 
of  structure  in  the  x-coordlnate 
direction 

Uniform  over  each  region  of  the 
surface  model  (pages  9,10) 

0PT5:  Reserved  for  later  use 


21-24 


=0000 


yoooo 


25-28 


=0000 


=0001 


29-32 


=0000 


OPT  6: 

Surface  geometric  data  (page  8) 
and  velocity  boundary  condi- 
tion data  are  read  in  from  cards 
for  each  station 

Automatic  generation  of  sur- 
face geometric  data 

0PT7:  Selector  for  near-field 

point  data 

Coordinates  (x-y-z)  of  near-  or 
intermediate-field  points  are 
to  be  read  in  from  cards 

Automatic  generation  of  near- 
or  intermediate-field  point 
coordinates 

OPT  8: 

Surface  pressures  and  near-  or 
intermediate-field  pressures 
are  to  be  calculated 
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0001 

Only  surface  pressures  are  to 
be  calculated 

0002 

Surface  pressures  and  far-field 
pressures  are  to  be  calculated 

33-36 

OPT9:  Selector  for  Induced  mass 
or  self-radiation  impedance  cal- 
culat ion 

0000 

Neither  quantity  is  computed 

0001 

Induced  mass  is  calculated 

0002 

Self-radiation  impedance  is  cal- 
culated 

37-40 

OPTIO:  Selector  for  far-field 

point  data 

0000 

Coordinates  (it;,*)  of  far-field 
points  (see  Figure  13)  are  to  be 
read  from  cards 

0001 

Automatic  generation  of  far-field 
point  coordinates 

SYMMETRY  OPTIONS" 

CARD 

Columns 

Contents 

Description 

1-  4 

0002 

yz-plane  (see  Figure  12)  symmetry 

0003 

yz-plane  antisymmetry 

CO 

lA 

0004 

xz-plane  symmetry 

0005 

xz-plane  antisymmetry 

9-12 

0001 

Longitudinal  (xy-plane)  symmetry 

0002 

Longitudinal  antisymmetry 

13-16 

0000 

No  rotational  (z-axls)  symmetry 

0001 

Rotational  symmetry 

'NUMBER  OF  REGIONS 

" CARD 

Co lumns 

Contents 

Description 

1-  4 

Total  number  of  surface  regions 
(pages  9,10,15)  used  in  modeling 
that  portion  of  the  structural 
surface  which  is  nonredundant 
with  respect  to  the  xy-plane 

• "REGION:  EXTENT  AND  MODELING" 

see  Figure  5) 


Columns 


17-2^* 


25-32 


Contents 


33-40 
41-48 
49-5  2 
53-56 
57-64 


=-1.0 

• "REGION:  EXTENT  AND  MODELING" 

surface;  see  Figure  10) 


CARD  (piecewise  conical  shell  surface; 


Description 

Radii  bounding  the  surface  region 


Angles  specifying  the  circum- 
ferential extent  of  the  sur- 
face region;  except  for  regions 
capping  the  ends  of  a body 
(9j  = 0 degrees,  6^  = 

360  degrees),  the  extent  is 
that  portion  of  a region  which 
is  nonredundant  with  respect 
to  principal  coordinate 
symmetry  planes 

Define  region  extent  on  the 

z-axis 

Control  automatic  modeling  of 

region  into  stations  (see  page  13) 

Parameter  designating  direction 
of  outward  normals  to  regions 
which  are  vertical  to  the  z- 
axis;  for  example,  rings  and 
disks  (Figure  6) 

Normal  points  in  z"*”  direction 

Normal  points  in  z direction 

CARD  (prolate  splieroid  and  spherical  shell 


Columns 


17-24 

25-32 


33-40 

41-48 


Contents 


Dcscript ion 


Blank 
Blan  k 


Same  as  given  for  these  parameters 
with  the  piecewise  conical  shell 
surface  case  above 
Define  extent  of  the  region 
in  colat 1 tude 
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% 


49-52 

"1 

Control  automatic  modeling  of 
region  into  stations  (see  pages 

53-56 

m| 

16,17,  Figure  11) 

"NORMAL  SURFACE  VELOCITY"  CARD 

Columns 

Contents 

Description 

1-  8 

V (Real) 

Real  part  of  the  normal  surface  in- 

n 

fluid  velocity  which  is  to  be 
assigned  each  station  (Figures  8 
and  11)  in  a region 

11-18 

V (Imag. ) 
n 

Imaginary  part  of  the  above 
velocity 

"PARii-METER  CARD  'A'  FOR  ASSUMED 

VELOCITY  DISTRIBUTION"  (see  page  20) 

Columns 

Contents 

Description 

1-  6 

Vp(Real) 

Real  part  of  constant  velocity 
amplitude  for  the  vibration  mode 

7-12 

v^dmag.  ) 

Imaginary  part  of  the  constant 
velocity  amplitude 

13-16 

M 

Integer  constant  associated  with 
the  vibration  mode 

=0 

Longitudinal  and  circumferentially 
symmetric  modes 

= 1 

Flexural  modes 

>1 

Lobar  modes 

* "PARAMETER 

CARD  'B'  FOR  ASSUMED  VELOCITY  DISTRIBUTION"  (see  page  20) 

Columns 

Contents 

Descript  ion 

1-  6| 

Rea  1 part  of  iKz  ) . 

1 

7-12  j 

Imaginary  part  of  il(z)^ 
where  ')'(z)^  designates  a 

13-18  1 

1 

discretized  point  of  the 

19-24  J 

1 

function  ^ilz)  assigned  to 
the  ith  circumferential  band 

■ 

of  the  surface  model  (for  ex- 

• 

• 

ample,  Figure  7) 

61-66 

1 

The  ordering  of  bands  is  from 

67-72 

I 

the  left  (z  ) end  of  the  body 

to  the  right  (z"*^)  end 


I 


i 


r' 


i 


K’! 


'V' 


1 1 


0 

1 ) 


FI 
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• "ELLIPSE  SEMI-AXIS  SPECIFICATION"  CARD 


Columns 
1-  8 
9-16 


Contents 

AA 

BB 


Description 
Semlmajor  axis 
Semlminor  axis 


• "SPHERE  RADIUS  SPECIFICATION"  CARD 
1-  8 


-•X  . 

1-  8 

9-16 

’ll 

17-24 

25-32 

- 

33-40 

ti’ 

41-48 

1 

49-56 

57-64 

V 

65-72 

1-  8 Radius  of  a spherical  surface 

• "SURFACE  GEOMETRIC  DATA"  CARD  (see  page  8) 

1-  8 Station  number 

9-16  Station  area 

17-24  \ X-,  y-,  z-coordlnates  of 

25-32  > station  base  point 

33-40  ) 

41-48  n^  X-,  y-,  z-coordinates  of 

49-56  n \ a unit  outward  (from  sur- 

y 

57-64  n face)  pointing  normal  vector 

z ' 

at  station  base  point 

65-72  Local  curvature  of  station 

surface  at  the  base  point 

• "IN-VACUO  SURFACE  MOBILITY  COEFFICIENTS"  CARD  (see  pages  5,8,9,19) 

Card  Columns  Contents  Description 

l-lb)  Real  part  of  q 


17-32 

33-48 

49-64 


Imaginary  part  of 

Real  part  of  q..  , 

-l- » j 

Imaginary  part  of  q.. , 

^ > J 

In  the  most  general  sense,  the 
mobility  coefficient  q^^ 

represents  the  normal  velocity 
response  at  station  1 due  to  a 
normal  unit  force  at  frequency 
u acting  at  station  j 

In  applications  Involving 
symmetric  surfaces  and 
symmetric  vibrational  load- 
ings, the  data  input  q^  may 
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represent  the  mutual  influence 
of  groups  of  stations  on  one 
another;  for  example,  In  the 
case  of  axially  symmetric 
surfaces  with  ring  loads, 


each  would  be  the 


velocity  response  of  band  (see 
Figure  7)  i due  to  a unit 
ring  load  at  band  j 


• "IN-VACUO  NORMAL  SURFACE  VELOCITIES"  CARD  (pages  5,19) 

For  fluid-structure  interaction  applications,  the  surface  normal 
velocities  for  the  set  of  stations  representing  the  nonredundant  portion 
of  the  surface  model  are  read  into  the  program  from  a tape  or  a series 
of  cards.  The  data  format  for  each  card  image  or  card  is: 


Columns  Contents  Description 

1-8  U (Real)  Real  part  of  In-vacuo  velocity  at 

n 

station  i 

11-18  U (Imag.)  Imaginary  part  of  in-vacuo 

n 

velocity  at  station  i 

The  ordering  of  the  velocity  input  is  according  to  the  station  numbering 
shown  in  Figure  9. 


Contents 


U (Real) 
n 


11-18 


U (Imag.) 
n 


• "AUTOMATIC  FIELD  POINT  GENERATION"  CARD  (pages  22,23) 


Columns 

Contents 

Description 

1-  3 

NFLAT 

Number  of  colatitudes  at  which 
field  points  are  to  be  generated 

4-  6 

NFLNC 

Number  of  longitude  points  to 
be  generated  at  each  colatitude 
specified  above 

10-17 

LATl 

Initial  colatitude 

18-25 

LAT2 

Final  colatitude 

26-33 

LONGl 

Initial  longitude 

34-41 

LONG  2 

Final  longitude 

42-49 

RADIUS 

Radius  of  spherical  surface  on 
which  field  points  are  to  be 

generated 
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• "ARBITRARY  FIELD  POINT  SPECIFICATION"  CARD  (pages  22,23) 

Columns  Contents  Description 


1-  8 
9-16 
17-24 


X-,  y-,  z-coordinates  of 
field  point 


• "AUTOMATIC  FAR-FIELD  POINT  GENERATION"  CARD  (pages  22,23) 


1-  3 


NFFLAT 


Number  of  colatitudes,  be- 
ginning with  0 degrees,  at 
which  far-field  points  are  to 
be  generated 


4-  6 


NFFLNG 


Number  of  longitude  points  to 
be  generated  at  each  colatitude 
(excepting  0 and  180  degrees) 
specified  above 


10-17 


LATLIM 


Limit  to  the  range  of  colati- 
tudes at  which  points  are  to 
be  generated 


18-25 


LONGLIM 


Limit  to  the  range  of  longi- 
tudes at  which  points  are  to  be 
generated  at  each  colatitude 
excepting  0 and  180  degrees 


• "ARBITRARY  FAR-FIELD  POINT  SPECIFICATION"  CARD  (page  22) 

Columns  Contents  Description 


1-  8 
9-16 


It, 


FF 


6, 


FF 


Colatitude  of  far-field  point 
Longitude  of  far-field  point 


The  precise  ordering  of  the  various  types  of  data  described  here  are 
indicated  by  the  input  forms  in  the  appendix.  At  present  there  are  two 
basic  orderings  depending  on  application. 

For  cases  involving  an  assumed  In-fluld  velocity  boundary  condition, 
the  data  on  form  1 is  followed  by  the  surface  model  and  velocity  data 
from  one  of  the  pairs  of  forms,  2,3;  4,5;  6,8;  7,8;  9,10;  or  form  11,  and 
then  data  for  the  near-  and/or  far-field  calculation. 
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For  cases  where  the  structure-fluid  interaction  effect  is  to  be  I 

calculated  and  the  In-vacuo  surface  velocity  and  mobility  data  are  on  I 

punched  cards,  these  data  are  placed  between  the  surface  model  data  | 

from  the  above-mentioned  pairs  of  forms  and  the  data  for  field  points.  | 

If  the  surface  velocity  and  mobility  data  are  stored  on  Internal  files  | 

(disk),  the  surface  model  data  are  followed  directly  by  the  field  data.  | 

EXAMPLE  CALCULATIONS 

CALCULATION  OF  FAR -FI  ELI)  SOUND  PRESSURES 

★ 

The  first  problem  is  taken  from  an  early  test  series  in  which 
XWAVE  was  used  to  calculate  far-field  sound  pressures  from  a radially 
vibrating  cylinder  with  fixed  end  caps.  The  calculation  is  performed 
for  the  case  kA  = 1,  kB/2  = 2,  where  we  select  k = 10  and  the  cylinder 
dimensions  A and  B as  shown  in  Figure  14. 


1 

I 

t 

I 


X 


Figure  14  - Geometry  of  Finite  Cylinder 


Since  there  are  three  principal  coordinate  planes  of  symmetry 
for  the  surface  geometry  and  velocity  boundary  condition,  only  one-eighth 
of  the  surface  is  specified  as  a basis  for  generation  of  the  total  surface 
acoustic  model.  The  specified  portion  of  surface  is  subdivided  into  two 
regions;  see  Figure  15. 


•k 

See  the  single  asterisk  footnote  on  page  3. 
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Figure  15  - Portion  of  Cylindrical  Surface  Subdivided  into  Regions 


The  region  modeling  data  used  as  input  for  XWAVE  specifies  that 
Region  1 shall  be  covered  by  a single  band  of  area  with  15  circumferential  j 

stations  and  Region  2 by  two  bands  each  with  15  stations  around.  This 
gives  A3  stations  covering  the  surface  in  Figure  15  and  through  their 
reflected  images  in  the  symmetry  planes,  a covering  of  the  whole  sirface 
by  8 X 45  = 360  effective  acoustic  elements. 

Because  the  surface  geometry  and  velocity  boundary  condition  are 
axially  symmetric,  the  surface  pressure  is  constant  for  each  band  and 
hence  calculated  for  only  a single  station  in  each  band.  The  station 
is  the  first  one  in  each  band  (see  pages  14,15,  and  19  (top))  as  in- 
dicated in  the  XWAVE  printout  for  surface  boundary  condition  and  pres- 
sures. Thus  for  this  problem,  we  are  required  to  solve  only  a system 
of  three  equations  for  the  surface  pressure.  Dimension  parameter  D2 
(see  page  24)  specifies  that  the  entire  matrix  will  be  in-core  during 
coefficient  generation  and  solution  time. 

We  request  that  the  iterative  process  terminate  when  either  a con- 
vergence criterion  0.0001  (pages  21,  22)  or  a limit  of  30  iterations  is 
met . 

J 
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ALL  APPLICATIONS 


DATA  INPUT  FORM  (1) 


PIECEWISE  CONICAL  SHELL  SURFACE 


DATA  INPUT  FORM  (2) 


FAR  FIELD  PRESSURE  CALCULATION 


NPUT  FORM 


ELAPSEO  TIME  AT  ENTRY  INTO  XWAtfE  = SECONDS 


flHllt  CTLIN0E9  • RADIUS.  A*0.1.  lEN6Ih.  B«0 


CDLArirUDC  « oc« 


NormallziriR  the  magnitudes  of  the  computed  far-fleld  pressures 
to  the  largest  magnitude  yields  the  curve  indicated  by  x in  Figure  16 
(reproduced  from  Figure  9 of  CMD-2ii-71).  The  solid  reference  curve 
shows  consistency  of  XWAVE  results,  using  increasingly  refined  surface 
models,  with  published  results  from  another  computing  method.^ 

CALCULATION  OF  SELF-RADIATION  IMPEDANCE 

The  second  example  is  an  application  of  XWAVE  to  calculate  the 
self-radiation  impedance  (see  the  section  entitled  OPTION  FOR  SELF- 
RADIATION IMPEDANCE  CAl,CLn.ATION)  of  a vibrating  ring  in  a spheroidal 
2 

baffle.  Figure  17  recalls  the  geometry  for  this  problem. 

The  ring  is  assumed  to  vibrate  radially  with  unit  velocity.  The 
velocity  over  the  baffle  is  of  course  zero.  In  the  original  investi- 
gation, calculations  were  made  for  a range  of  values  of  the  acoustic 
shape  parameter  h = kd/2.  For  the  illustrative  calculation  given  here, 
we  choose  a value  of  3.0  for  h;  this  corresponds  to  the  wave  number 
k = 1. 374773. 

As  in  the  first  example,  there  are  three  principal  coordinate 
planes  of  symmetry  for  the  surface  geometry  and  boundary  condition.  The 
specification  of  the  poition  of  surface  that  is  to  be  divided  into 
regions  and  stations  proceeds  slightly  differently  in  this  case  because 
the  ends  of  the  spheroid  are  modeled  by  regions  consisting  of  single 
stations  centered  about  the  respective  poles;  see  Figure  18.  The  data 
specify  the  end  cap  region  extending  from  e = 0 to  fl  - 360  degrees 
rather  than  to  9 = 90  degrees  as  with  Regions  2 and  3 since  the  only 
imago  data  obtainable  from  the  basic  geometry  (for  example,  coordinates 
of  base  point  and  unit  surface  normal,  etc.)  of  the  single  station 
within  the  region  arc  for  the  opposite  end  through  longitudinal  re- 
f lection . 


Schenck,  H.A.,  "Improved  Integral  Formulation  for  Acoustic  Radiation 
Problems,"  Journal  of  the  Acoustical  Society  of  America,  Vol.  44,  No.  1, 
pp.  41-58  (Jul  1968). 
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Figure  18  - Surface  Regions  for  Spheroidal  Baffle  with  Ring 


I 


The  data  for  modeling  Region  2 indicate  a covering  by  eight  bands, 
each  with  eight  circumferential  stations,  and  a covering  of  Re,  3 
(which  extends  one-half  the  ring  width)  by  three  bands,  each  with  22 
stations. 

The  effective  covering  of  the  entire  surface  is  then  8 x (8x8+3 
X 22J  + 2 (end  caps)  = 1042  acoustic  elements. 

Taking  account  of  radial  symmetry  of  the  surface  geometry  and 
boundary  condition,  the  program  generates  the  coefficients  for  a system  of 
12  equations  for  surface  pressure  (12  is  the  total  number  of  bands  includ- 
ing the  end  cap).  The  iterative  solution  is  accomplished  with  dimension 
parameter  D2  specifying  no  more  than  two  rows  of  the  coefficient  matrix 
in-core  at  each  stage  of  the  process. 

Option  9 (page  33)  selects  the  radiation  Impedance  calculation. 


/ 
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ALL  APPLICATIONS 


PLANE  PLANE  LONG.  RADIAL 


PROLATE  SPHEROIDAL  SHELL  SURFACE 


yiBRtTlHG  ^<IRG  IX  SI>HERalD>L  BAFFLE 
SElF-B»OI«TION  IXPEOANCE  OF  RING  FOB  H»3 


9rr 


9.S?1U-Q? 


NUMERICAL  CALCULATION  OF  INDUCED  MASS 


The  third  and  final  problem  illustrates  a numerical  calculation  of 
Induced  mass  according  to  the  method  and  formulation  discussed  In  the 
section  entitled  OPTION  FOR  INDUCED  MASS  CALCULATION.  The  problem  posed 

here  is  based  on  an  earlier  numerical  Investigation  of  the  Induced  mass 

* 

effect  for  cylinders  In  rigid  body  motion. 

We  consider  a cylinder  of  length  4.11  cm,  radius  1.27  cm,  and 
oriented  with  respect  to  coordinates  as  the  cylinder  in  Figure  14.  The 
cylinder  vibrates  perpendicular  to  its  longitudinal  axis  with  assumed 
velocity  magnitude  unity.  This  motion  generates  a normal  surface 
velocity  pattern  which  is  symmetric  with  respect  to  the  xz-  and  xy-planes 
and  antisymmetric  with  respect  to  the  yz-plane.  We  therefore  specify, 
as  in  the  first  example  problem,  only  an  octant  of  the  surface  as  a 
basis  for  generation  of  the  acoustic  model.  Figure  15. 

For  the  present  calculation,  the  data  specify  a modeling  of  the 
end  cap.  Region  1,  by  six  bands  with  nine  stations  per  band  and  of  the 
lateral  surface.  Region  2,  by  12  bands  each  with  nine  stations.  This 
gives  an  effective  covering  of  the  entire  surface  by  8 x (6  x 9 + 12  x 9) 
= 1296  elements.  Thus  for  our  purpose,  we  are  using  here  a considerably 
less  refined  surface  model  than  the  3456-element  model  used  to  obtain  the 
previous  results  (CMD-41-71). 

The  input  value  k.  = 0.00001  is  an  approximation  for  the  condition 
of  fluid  incompressibility.  The  fluid  density  is  unity. 

Option  4 (page  32)  is  selected  to  obtain  automatic  generation  of  the 
normal  surface  velocity  at  each  station  from  the  velocity  magnitude 
specified  for  each  region  on  the  "NORMAL  SURFACE  VELOCITY"  data  cards. 

In  the  case  of  rigid  body  motion  as  considered  here,  the  velocity  data 
on  these  cards  do  not  signify  normal  surface  velocities  but  rather 
velocities  in  the  direction  of  rigid  body  motion,  from  which  normal 
components  will  be  automatically  generated. 

The  induced  mass  (in  grams,  since  centimeter-gram-second  system  was 
used  in  data)  obtained  here  is  within  2.7  percent  of  the  results  ob- 
tained with  the  finer  model. 


See  the  double  asterisk,  footnote  on  page  3. 
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DATA  INPUT  FORM  (1) 


PIECEWISE  CONICAL  SHELL  SURFACE 


DATA  INPUT  FORM  (2) 


OIMCNSIOMS  fOR  «99«rS 


l.tOb«-ai  -l.?9eir»10  ?.9<.00t*00  i.79^9^-a^  9.9tl9C-01  «.7l9»£-0? 


1915E-I1 


^^87E♦9I>  J.ESEOf-Ol  -9. %t Jtf-0 1 ?.5li05E»01!  !.7<^5<»E-0^ 


enter  SUBROJTlNf  EOR  IT£9ATtVE  SOLUTION  fOR  SUREtCI  PRESSURE 


S0T7P-18 


5l^E-lt 


-7.9W35C-0b 


S135E-13  5hEI>« -OS  -?.CSS9£-U  -5.qa??t-0fc  -^.S3^le-l3  -7.3I77E-06  -^.93  3 3£-l3  -O.SESkE-tO 
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DATA  INPUT  FORM  (1) 


PIECEWISE  CONICAL  SHELL  SURFACE 
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(CONTINUATION  SHEET) 


PROLATE  SPHEROIDAL  SHELL  SURFACE 


SPHERICAL  SHELL  SURFACE 


PROLATE  SPHEROIDAL  SHELL  SURFACE  OR  SPHERICAL  SHELL  SURFACE 

(CONTINUATION  SHEET) 


PROLATE  SPHEROIDAL  SURFACE  WITH  ASSUMED  VELOCITY  DISTRIBUTION  (SEE  PAGE  20) 


♦ CONTINUATION  CARDS  AS  NEEDED  DATA  INPUT  FORM  (9) 


TITLE  DATE 


ARBITRARY  SHELL  SURFACE 


DATA  INPUT  FORM  (111 


STRUCTURE  FLUID  INTERACTION  APPLICATIONS 


INPUT  FORM  (12) 


STRUCTURE  FLUID  INTERACTION  APPLICATIONS 


INPUT  FORM  (13) 


DATA  INPUT  FORM  (14) 


DATA  INPUT  FORM 


INITIAL  DISTRIBUTION 


Copies 

1 CNO  OP  09  8T 

1 CHONR /Library 

I CHONR/N.  Basdekas 

Code  439 

1 CHONR  439/N.  Per rone 

1 DNL 

1 USNA/Dept  Math 

1 USNA/Lib 

1 NAVPGSCOL/Lib 

1 NROTC  & NAVADMINU,  MIT 

1 NAVWARCOL 

1 NAVSEA/CDR.  Cox 

1 NAVSEA  037/C.  Taylor 

1 NLONLAB  NUSC/A.  Carlson 

1 NAVSHIPYD  PUGET /Lib 

1 NAVSHIPYD  CHASN/Llb 

1 NAVSHIPYD  MARE /Lib 

1 NAVSHIPYD  NORVA/Lib 

1 NAVSHIPYD  PEARL /Lib 

1 NAVSHIPYD  PHILA  240 

1 NAVSHIPYD  PTSMH/Llb 

12  DDC 


CENTER  DISTRIBUTION 


Copies 

Code 

1 

1720 

M.A.  Krenzke 

1 

1720.3 

R.F.  Jones,  Jr 

1 

1720.3 

T.E.  Reynolds 

1 

1800 

G.H.  Gleissner 

1 

1805 

E.H.  Cuthill 

1 

1840 

H.J.  Lugt 

1 

1843 

J.W.  Schot 

1 

1843 

C.W.  Dawson 

1 

1844 

S.K.  Dhir 

1 

1844 

G.C.  Everstine 

20 

1844 

F.M.  Henderson 

1 

1850 

T.  Corin 

1 

19 

M.  Sevik 

1 

1901 

M.  Strasberg 

1 

1902 

G.  Maidanik 

1 

1903 

G.  Chertock 

1 

1960 

D.  Feit 

1 

1962 

A.K.  Kukk 

1 

1962 

A.  Zaloumis 

1 

1965 

J.  Brooks 

1 

1965 

H.B.  Ali 

1 

1965 

Y.N.  Liu 

1 

2740 

L.J.  Argiro 

1 

2740 

Y . F.  Wang 

30 

5214.1 

Reports 

Distribution 

1 

522.1 

Library  (C) 

1 

522.2 

Library  (A) 
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tt)  OtWSRW  REPORTS,  A FORMAL  SERIES  FURLWHINO  iRfORMATlOH  OF 
FBRMANENT  TeORlKAL  VALUE,  OE^  A SERIAL  REPORT  NUmER. 

12)  OfP>lRtl(l»ITAL  REFbR|».  A SEMIFORMAL  SERIES,  RECORDING  INFORMA 
TION  Of  A I^LIMINARY  OR  TEMPORARY  NATURE.  OR  OF  LIMITS  INTEREST  OR 
significance.  CARRYING  A DEPARTMENTAL  ALPHANUMERIC  IDENTIFICATION. 

(3)  TBCMNiaa  MEMORANDA.  AN  INFORMAL  SERIES,  USUALLY  INTERNAL 
WORKING  PAPERS  OR  DIRECT  REPORTS  TO  SPONSORS.  NUMSEREO  P3  TM  SERIES 
R0ORTS;  NOT  FOR  GENERAL  DISTRIBUTION. 


